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Abstract 

A brane-antibrane model for the entropy of neutral black branes is developed, fol- 
lowing on from the work of Danielsson, Guijosa and Kruczenski [1]. The model involves 
equal numbers of Dp-branes and anti-Dp-branes, and arbitrary angular momenta, and 
covers the cases p = 0,1,2,3,4. The thermodynamic entropy is reproduced by the 
strongly coupled field theory, up to a power of two. The strong-coupling physics of the 
p = case is further developed numerically, using techniques of Kabat, Lifschytz et al. 
[2,3], in the context of a toy model containing the tachyon and the bosonic degrees of 
freedom of the DO-brane and anti-DO-brane quantum mechanics. Preliminary numeri- 
cal results show that strong-coupling finite-temperature stabilization of the tachyon is 
possible, in this context. 
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1 Introduction 



The drive to explain the thermodynamic entropy of black holes and black branes via the 
statistical mechanics of microscopic degrees of freedom has been a central preoccupation 
of string theorists and other gravitational theorists for at least two and a half decades. 
String theory has provided significant progress on this front in the last eight years. The 
celebrated 1996 success of Strominger and Vafa in computing the entropy of particular 
D = 5 supersymmetric (BPS) black holes occurring in low-energy string theory, by 
using a conformal field theory appropriate to the microscopic physics of strings and 
D-branes, even raised the profile of string theory itself in a significant way. Further 
successes followed; the thermodynamic entropy of BPS and near-BPS black holes in 
various dimensions was reproduced successfully, in some cases even beyond leading 
order in macroscopic quantum numbers such as mass, charge, and angular momenta. 
A microscopic string theoretic accounting of the entropy of neutral black holes and 
black branes remains more elusive, however. One reason is that neutral black holes 
and branes are literally as far from BPS as possible, meaning that supersymmetry may 
provide no help at all in the endeavour. 

BFSS matrix theory |3| is a conjectured relationship between quantum D = 11 
M theory in light-front frame and the supersymmetric quantum mechanics of a large 
number N of DO-branes. This theory was recruited, e.g. in j^j, to help explain the 
entropy of near-BPS black holes in dimensions D = 10 — p, and also in attempts to 
explain the entropy of the neutral cases as well. Reviews of this story include [0] (see 
references therein). Another interesting proposal was in [7j. Here, we simply record 
one major common feature of these various matrix theory approaches to counting the 
entropy of neutral black holes, which is relevant to our work here. Namely, the need 
for an infinite boost in the 11th dimension x^^ in order to relate the neutral black holes 
of interest to nearly-BPS systems whose entropy can be computed microscopically 
using D-brane field theory. In particular, this infinite boost in the 11th dimension 
eliminates anti-DO-branes from the picture. The approach that we will take here will 
be different. Other interesting approaches to neutral black hole entropy from quite 
different perspectives include jHj and references therein. 

Gravity/gauge correspondences are specific relationships between open string and 
closed string degrees of freedom. For the case of sixteen supercharges jH], they relate 
the degrees of freedom of supersymmetric Yang-Mills theories to closed string theory 
on near-horizon Dp-brane spacetimes. Subsequent string theoretic developments of 
interest for our work here include tachyon condensation in brane-antibrane systems. 
The paper which sparked our specific interest was J,, , in which an interesting proposal 
was made for a microscopic explanation of the entropy of neutral black D3, M2 and M5 
(non-dilatonic) branes in terms of a system of strongly coupled branes and antibranes. 

In section 2, we review features of earlier analytic work, upon which we build here. 
Section 3 contains the bulk of our analytic observations. We study the entropy of 
neutral black Dp-brane spacetimes by using a microscopic model with an equal number 
of strongly coupled Dp-branes and anti-Dp-branes, with arbitrary angular momenta 
turned on. Our main result is that the supergravity entropy density can be reproduced 
by strongly coupled field theory, up to a power of two. Section 4 contains an exposition 
of preliminary numerical work on the p = system. Our primary goal in the numerical 
part of our work is to try to test some assumptions of the brane-antibrane model, 
by doing direct strong-coupling simulation of the p = system, adapting methods of 
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[Slini. We use a bosonic toy model to do this analysis, and find preliminary evidence 
that tachyon stabilization may indeed occur in the fashion expected from the work of 

We are not yet able, however, to be conclusive about the supersymmetric case. 

While this paper was in preparation, the work of j2J appeared, which has some 
overlap with section 3, for the case p = 3 with angular momentum. 

2 Analytic ingredients from prior work 

2.1 Gravity/gauge correspondences for general p 

For a system of Dp-branes, it is possible^ to take a clean low-energy limit such that 
the open strings with endpoints on the Dp-branes decouple from the closed strings in the 
bulk. This observation led, of course, to the celebrated AdS/CFT correspondence and 
the non-conformal open/closed-string correspondences of Itzhaki, Maldacena, Sonnen- 
schein and Yankielowicz The theory on the branes is supersymmetric Yang-Mills 
theory (SYM) with sixteen supercharges. The 't Hooft coupling for the SYM theory is 

4MA^^(2^rvr'iv, (1) 

where gs is the string coupling and £s = x/c? is the string length. 

These open/closed string correspondences tell us that we can hope to understand 
the gravitational fields of Dp-branes via the dual SYM theory. For p ^ 3 the fields 
describing the D = 10 geometry vary radially in the spacetime, and on the SYM 
side the coupling is dimensionful, yielding breakdown of the weak-coupling description 
either in the UV or the JR. In other words, there are limits to the validity of both 
descriptions. A dimensionless control variable is given by g'^g = ^ym-^^^"^' 
supergravity side, U = r/i^ is the radial isotropic coordinate in energy units. The 
requirements [5] for the supergravity geometry in = 10 to remain valid are 

1 < {g^^NUP-^) < . (2) 

At the left-hand end, a' corrections become important and the SYM theory takes over 
from the D = 10 geometry as the weakly-coupled description; at the right-hand end, 
strong coupling (dilaton) ensues and it is necessary (for large- A^) to turn to the S-dual 
supergravity geometry. 

It was earlier noted in 1(7 that thermodynamic properties of near-BPS Dp-brane 
(and M-brane) geometries could be written in a way that is suggestive of a field theory 
interpretation. The main focus of that paper was the non-dilatonic branes, for which 
the energy density above extremality Am and entropy density s were written in a 
way reminiscent of gases of weakly interacting massless particles in d = p + 1. For 
the dilatonic cases, however, the expressions do not yield recognizable weak-coupling 
results. The physical interpretation of this fact is that the strongly coupled physics 
of the p ^ 3 SYM theories gives rise to nontrivial dependence of the entropy and the 
energy above extremality on the temperature, 

2(7-p) (9-p) 

Am(r) oc T , s(r) oc r(5-f) , (3) 
For p < 4, which are the cases on which we will concentrate 
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where T is the Hawking temperature of the geometry. 

In the p = case, checking aspects of the open/closed-string correspondence con- 
jectures is potentially feasible. The reason is of course that the SYM theory in this 
case is actually matrix quantum mechanics - albeit with sixteen supercharges. The nu- 
merical investigations of [21 Ej aimed to check the p = correspondence explicitly, by 
finding the entropy of the system of strongly coupled N DO-branes and comparing it to 
the entropy of the near extremal DO-brane supergravity background. Results obtained 
showed that the 9/5 power in the expression for the entropy in © was indeed approx- 
imately reproduced in the numerical approach. This striking result, in combination 
with the work given a lightning review in the next subsection, provided the essential 
motivation for our work. 

2.2 Brane and antibranes at finite temperature and the 
entropy of black branes 

Motivated by an observation of Horowitz, Maldacena, and Strominger Daniels- 
son, Guijosa and Kruczenski ^ argued that, starting with a brane-antibrane system 
at zero temperature, turning on a finite temperature could lead to reappearance of 
open string modes. This physics was argued to ensue at strong open-string coupling 
QsN but weak closed-string coupling gg, as would be appropriate to a model with va- 
lidity in the black brane regime. Note that, to have this reappearance of the open 
string modes at temperatures below the Hagedorn temperature, going to the regime of 
strong open-string coupling was crucial. In alternative language, the tachyonic mode 
becomes stabilized by finite-temperature and strong coupling physics. Based on these 
observations, the authors of T formulated a microscopic brane-antibrane model for 
the entropy of non-dilatonic black branes: the D3, M2 and M5 cases. 

The model has two SYM theories on world volumes, one for the branes and the 
other for the anti-branes. The tachyon, argued to be stabilized by finite-temperature 
and strong coupling effects, does not contribute to the total energy of the black hole, 
because it sits almost on top of the tachyon potential. It is also argued that the tachyon 
does not contribute to the total entropy, because the tachyon gets a large thermal mass. 

The black holes of interest to us are neutral. Since the black 3-brane in D = 10 
(or, equivalently, the black hole in D = 7) is neutral, it is taken to be modelled by 
the gas on a number of D3-branes, the gas on an identical number D3-branes, 
and the tachyon dynamics. Since the branes and anti-branes appear in equal numbers, 
it makes sense to assume that the corresponding temperatures are equal. Using these 
assumptions, for the strongly coupled field theory side, 

mpT = 2Nt3 + a—N^T^, 

o 

SFT = a^N^T^ , (4) 



where T3 is the D3-brane tension and a is a constant known from AdS/CFT to be 6 
(not 8) when the SYM theory in question is strongly coupled. Note that the quantities 
of thermodynamic interest written here are specific, i.e. mass and entropy densities. 
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Rearranging to eliminate T gives 



.PT = a(vr/6)iV ^-^-.^^J . (5) 

It is important to emphasize that this analysis is done in the microcanonical ensem- 
ble, where the total mass and charge of the black hole, as well as any angular momenta 
present, are kept fixed. Working in the microcanonical ensemble prevents thermal 
fluctuations from creating an infinite number of D3-D3 pairs. Working in canonical 
picture, by contrast, would allow the system to catastrophically create an infinite num- 
ber of pairs, by extracting an arbitrary amount of energy from the reservoir, as it is a 
thermodynamically favourable state. In other words, one indefinitely spends energy to 
make more entropy, which makes the thermal ensemble destabilized. 

The essential innovation in ^ is that the total number of branes is regarded as 
an independent variable, which is determined thermodynamically. The next step in 
the analysis, then, is to maximize the entropy with respect to A^. The value of N so 
obtained is related to the mass density as 

mpT = 5A^T3 . (6) 

To proceed further, it is necessary to decide what to do with the masses. One 
assumption is to take mpT = "^SG- Using this, and substituting back for the optimized 
value of leads to the final expression for the entropy from the strongly coupled field 
theory side, 

1 5 _1 _5 1 _ I 

SFT = a*243 45 4Tr4^nrn^^, (7) 
where k = \/^/t-^. One the other hand, entropy of a black 3-brane is 



. 9 _ 5 1 
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SSG = 2*5 iTT-i^/Km^Q. (8) 
Identifying msc = mpT and a = 6 gives 

ssG = ■ (9) 

This was the result of pQ. It shows that the entropy scaling from the brane-antibrane 
model is correct. It is worth noting here that this scaling agreement is nontrivial: it 
does not come about through simple dimensional analysis (and counting powers of N). 

The overall coefficient misses, by a factor close to unity. It is perhaps not surprising, 
though, that the agreement, while close, is not exact. Reasons for this may include the 
fact that there is no clean decoupling between open string and closed string degrees 
of freedom. These effects not taken into account in the brane-antibrane model may 
indeed play a role, but they cannot be major effects for quantities like the entropy, 
because the scaling comes out correct. 



3 Dp and D^^-branes with angular momenta 

In order to test the idea of the brane-antibrane model for neutral black hole entropy 
further, it is interesting to consider the cases other than p = 3, and to add angular 
momentum. This will be the main focus of the analytic part of our work. 
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It is convenient to begin by reviewing some salient properties of rotating black Dp- 
branes for various p. We are particularly interested in the p = story, since that is the 
case of the brane-antibrane model which we plan to test numerically, starting with the 
preliminary investigations of section 4. 



3.1 Supergravity side 

Black p-branes in D = 10 are of course equivalent to D = 10 — p dimensional black 
holes. We will not need the precise form of the supergravity fields; what is important 
for us here is the relationship between various physical parameters including tension, 
charge, angular momenta, and horizon radius. Using for example JHI, and converting 
to quantities which are specific (per unit volume), we have for the D = 10 black branes 



Th 



m 



(7 - p) levrG ° 
(7 - p) - 2k 
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4TTrH 
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dr. 



16ttG 



(10) 



where C is the boost rapidity parameter while rQ,li are the original parameters used 
in creating the solutions, while the Newton constant is given by IGnG = (27r)^(7g^^. 
We also use the shorthand ap = (7 — ^)^8 
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Values of ap are tabulated here for 



convenience 



p 
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57r^ 
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Stt 


27r 



The parameter rn is the horizon radius given by 



(11) 



It is straightforward to massage the expressions to write the entropy density s in terms 
of the energy density above extremality Am. To do that, it is convenient use (|llj) to 
obtain tq in terms of th and 



This is simple if we make the definitions 

A,. A. 



Then we have 



11(1 + ^? 



1 



(12) 



(13) 
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The energy density above extremality Am = m — q\s given by 



Am = — ^ 
levrG ° 



1 1 



:i + c) (7-p) 

We can also write the Hawking temperature as 



(14) 



th = v^' ' (15) 



47rro p{\i) 

where p{\i) is given by and we use the additional shorthand 



9 " \2 



(7-/.)fe(l + Af) 

In addition, the number of Dp-branes can be written as 



where 

1 



(18) 



" <7s(2^X+^ 

is the Dp-brane tension. 

Our main interest is to explain the entropy of neutral black branes m. D = 10, or 
equivalently, neutral black holes in D = 10 — p. We therefore record the expressions 
here that we aim to reproduce using the strongly coupled brane-antibrane field theory 
model. 

In the case of a neutral black hole, we have Am = mgc- Defining 



A 

(7-p)' 



^ - )^ > (19) 



we have for the energy density 

— = life--""'*' 

and for the Hawking temperature 



{7-p) kp{\) 

^^-^^7(aT- ^''^ 

Therefore, 

_«p A A7-p) fcp(Ai) 1 

while for the entropy 

« "P „t ^ ^ ( - P) hi^i) \^~'^ /9o^ 

'^^^^^^ - (7^i6^^^^^^ [r^Ti^T^i ) • ^^^^ 



Defining 

= 1 16^J 
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(24) 



we have that 



ssg("t-sg) 



27r 



(25) 



We may now ask how to unfurl the dependence of p{Xi) on msG cind the rotation 
parameters ji. The simplest way to proceed is to recognize that 



K _ ji 
2tt~ 1 



(26) 



We therefore have the equation 

n 



.i=l 



(7-p) 



1 

J- />-vi 



SG 



where 



Then, finally, we have 



ssG = 2-np{RsG (msG , G) , ji ) Rsg ("t-sg , G) . 



(27) 
(28) 

(29) 



We want to invert H27|) to unfurl the dependence of Aj on msG,G, ji. When no 
angular momenta are present, we have of course that p = 1. In general, however, the 
equation for Aj is a polynomial of degree (7 — p) in A, which is a quintic or worse for 
p <2. For now, we do not concern ourselves with whether we can actually solve for p; 
we just leave the dependence of p on (mgGi G,ji) implicit. We also note that for p = 3 
and one angular momentum parameter ji, we can actually solve to find 



and 



where 



ai(:p = 3) = ./2^^+yr+^ 



pip = 3) 



(30) 



(31) 



(32) 



3.2 Near-extremal: the field theory side 

In this subsection, we take the near-extremal limit of various supergravity formulae to 
tell us how the strongly coupled field theory quantities should behave on the branes 
and anti-branes, using IMSY duality. 

We should mention that in some respects our analysis here is quite similar to |lUj . 
There are two major differences, however. The first is that we are interested in the 
effect of turning on angular momenta ji. Our results must of course reduce to those 
of ^U] upon taking ji = 0. The second is that we plan to use the information about 



7 



the strongly coupled field theory, obtained via weakly coupled supergravity, using the 
techniques of p. 

Using the relation (|17|) , we can^ express tq in terms of Np 



7-P 



2(1 - C) NpTp- 



IQttG 



OLr 



(33) 



(Of course, this tq is appropriate to the near-extremal geometry, and is not the same 
as the ro at the end of the previous subsection, which refers to the horizon radius of 
the neutral spacetime.) Equivalently, 



2(1-0 = 



The energy density above extremality becomes 



Am^r 



(1 - ONpTp 



(9-p) 



Near extremality we have^ for the Hawking temperature 



(34) 



(35) 



(5-p) 

Th = [2(1-C)]^(^ 



± f 
p) - 



OLr, 



V levrGiVpTp 



(5-p) 



[2(1 - C)(27r)^]^(^ 



(27raj; 



{7-p)kp{Xi) 
4tt p{\i) 

(7-p) kp{\i) 
47r p{\i) 



(36) 



We can obtain the equation of state for the near-extremal black Dp-brane by eliminating 
Q in favour of Th , 



Ambranes(7H) = {-fp{27r f {27rap) (^-p)}x 

4tt p{\i) 



{7-p)kp{Xi) 



Th 



2(7-p) 
(5-p) 



Defining the abbreviation 



7p 



(9-p) 
2(7 -p) 

we have for the entropy density on the branes 



(37) 



(38) 



•S branes 



(7-p) 



(27r)^ (27rap) (s-f) \ x 



(p-3) 



47r p{\i) 



Th 



(5-p) 



(39) 



(7-p) fcp(Ai 

This equation is of central importance; it encodes the equation of state for the system. 



^The constant Cp familiar from the TASI-99 lectures of one of us is related to the quantities used here by 
Cp9stl~^ = 167rG/ap. 

•^For reasons discussed in etc., the cases p = 5,6 are more problematic to interpret. We therefore 
restrict ourselves from here on to the cases p < 4. 
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Now we come to the crucial step. Following the innovation in we actually take 
the lead for the behaviour of the strongly coupled field theory on the branes and anti- 
branes by using the near-extremal supergravity results. This is tantamount to using 
IMSY duality Here, for general p, our proposal to use IMSY duality for the p ^ 3 
case in the model of type ^ is a more nontrivial step than in the p = 3 case where 
the field theory behaved in a simple way. This assumption involves some nontrivial 
physics; like pQ, we are not taking into account the lack of a clean decoupling limit 
in the brane-antibrane model. Nonetheless, it is our point of view here that taking 
the near-extremal supergravity result seriously for the strongly coupled field theory on 
both the set of branes and the set of antibranes is exactly what we need. 

We are particularly interested to investigate this story for the DO-DO case. The 
reason is that in the d = + 1 case we have the hope of actually checking the above 
assumption explicitly, using a strong-coupling simulation. In particular, a significant 
motivation for our investigation of the non-conformal (p ^ 3) cases in the first place 
was the result of |21 E] in which the equation of state was approximately reproduced 
numerically. We will begin developing the numerical story for the DO-DO case in the 
next section, but for now we work out the analytics. 

For the open-string gas on (say) the set of Dp-branes, we have 

Sbranes = J-p{Xi)^9^ (Am)> , (40) 
Op 

where we have defined 

6p-^7?(2^)^ap^. (41) 

Now, for our model, recalling that we have the field theory on both the branes and 
antibranes, we have 

mpT = {2)NTp + (2)Am(rFT) • (42) 

BPS branes carry no macroscopic entropy, so we take the entropy of the strongly 
coupled field theory system representing the neutral black brane to be twice the entropy 
of the gas on each set of branes, sft(2ft) = (2)s(7ft)- We find 

SPT = {2f^piX^)^g^ - Nrpy^ . (43) 

We now need to know how p{Xi) depends on variables of interest. We again use the 
trick of the previous subsection to write Aj = 27rjj/s. Then 

1 

A.p(A.) = AJn(l + A?)| = ^^^(9YM)^J^[^-Nrp^ 



-7p 



ji 

Rft{N, mFTiffYM 



where 



RFT{N,mFT,gYM,rp) = {2)yVNg^^''^ V~(2)~ ~ ^'^f 



, 7p 



(44) 



(45) 
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Referring back to the entropy equation (|43)) , we have for general p and general angular 
momenta that 

SFT = 2TTRFTiN,mFT,gYM,Tp)p{RFT{N,mFT,gYM,Tp),ji). (46) 

At this stage, we may wonder whether the equations for the physical parameters 
Xi in terms of {ji, ruFT, N, qym, Tp) can actually be solved explicitly analytically. For 
some p, they can. For other cases, including the p = case of most interest to us, 
however, they cannot. For now, we will put this issue aside, and just proceed with the 
dependence of p on Rft and thereby on {N, itift, gym, Tp, ji) implicit. 



3.3 Comparing field theory and supergravity 

We now take the field theory result of the last subsection and ask what happens when 
we optimize with respect to N, the number of branes (and anti-branes). For simplicity, 
we first ask how this story works without angular momenta. Later we add angular 
momenta back in. 

Optimizing the entropy of the strongly coupled brane-antibrane theories w.r.t. N 
gives 

^= 9 n7'? \ mFT = {2NTp)26p. (47) 

2Tp[l + 2-fp) 

The energy density in the brane gases is then 

mpT - 2NTp = {2NTp)2-fp . (48) 

Substituting these expressions for N and mpT back into the expression for the entropy, 
and converting Dp-brane quantities into the Newton constant, we obtain 

271- c 

SFT = 2-> — G^f-^mp^T • (49) 
dp 

Two possible conclusions can be drawn from this. The first is that 

niFT = iTT-SG and sft = 2 ^C^-p) ssg (50) 
Alternatively, we can conclude that 

(9-p) 

Sft = Ssg and itift = 22(8-p) msc ■ (51) 

The second interpretation has an interesting conclusion. It says that the field theory 
energy, which is simply the energy on the Dp-branes plus the energy on the Dp-branes, 
is not simply the supergravity energy, but there is a (suitably negative) binding energy 

ImbindingI = msG (^2W^ - 1^ . (52) 

We find this interpretation the more attractive one. Binding energy can be expected 
in our model, because of the lack of a clean decoupling limit in our system between the 
open-string and closed-string modes. 

Let us now add back the angular momenta and see if it affects our exposition of 
the basic physics of the brane-antibrane model. 
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Recall that, in the micro canonical ensemble that we are studying, both mpT and 
ji are constants. This is why we have labelled i?FT(-^) a function of N, which we have 
to optimize following the model of 

Regardless of whether the equation for Aj can be solved explicitly, the entropy 
depends only on /9(i?FT)j as in (|l6|) . So we just proceed with p{Rp^) implicit. 

Our principle is to maximize the entropy as a function of A^. We have 

1 dspT 1 dRpT ^1 dp dRpT 



SFT dN i?FT dN p{RFT)dRFT ON 



1 I dp 

+ 



i?FT p(-Rft) 9-Kft 



Since we do not in general know the analytic dependence of p on i?FT) we need to find 
the derivative implicitly. We have from ()45|) 

< _ ^ 

^'pdR~ R^p dR- 

Using (|i6]) . and Aj/(27r) = ji/s, we have Aj = ji/{pR), which gives 



1 1 dp 

+ 



dRpT 



i?FT p{Rft) ORft 
Referring back to (@U), it is easy to find dXi/dR, 

[(7_p) + (5_p)A2]A^^ dX. _ (7-p) 



(55) 



(1 + A?)n,(l + A,2) 9i?SG i?^G''^'' ^^^^ 

Therefore the derivative is manifestly negative, as long as p < 4, and it is nonzero as 
long as at least one angular momentum parameter (jj) is turned on. Consequently, the 
term in square brackets in (|53() does not vanish, and it is nonsingular. Therefore, what- 
ever the behaviour of p{Rft), we are safe in concluding that the entropy is extremized 
by demanding that dRFT{N)/dN = 0. Optimizing Rft{N) we find that 

N= (57) 

which is exactly what we had for the non-rotating case. 

Substituting back this optimal value of into the field theory quantity Rft, we 
find that 

RFT{mFT,jiFT,9YM,Tp) = Rscim, jiSG, G) , (58) 

where Jift are the angular momenta in one copy of the strongly coupled field theory, 
if we make the identifications 

(9-p) 

mpT = msG22{8-p) , jiSG = jiFT ■ (59) 

Now, naively we would have expected the total angular momenta of the neutral su- 
pergravity solution to be split in half, shared equally between the strongly coupled 
brane and antibrane field theories. The fact that the angular momenta and the mass 
do not match precisely is not particularly surprising, however, because of the lack of 
a clean decoupling limit. In fact, closed-string modes (whose physics not included in 
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the brane-antibrane model) might carry both mass and angular momenta. We find it 
intriguing, though, that the renormalization factors are simply powers of two! 

Now, in order to get the entropy to match between the field theory side and the 
supergravity side here, we have the freedom to apply renormalization factors to both 
the mass and angular momenta. Alternatively, there is insufficient information to set 
these renormalizations factors unambiguously using our analysis. 

Let us then accept the mass and angular momenta renormalizations of (|59|) . We 
then find the remarkable fact that the functional form of the field theory entropy (|4H|) 
in terms of R is identical to the functional form of the entropy on the supergravity 
side (j29j) . With the renormalizations we see that the i?'s are the same on both 
sides. The strongly coupled field theory of branes and antibranes therefore reproduces 
the supergravity result, up to renormalization factors of two, i.e. 

/ (9-p) \ 

SFT ( "ipT = 22(8-p) mscJiFT 1 = SsG{msG,jiSG) ■ (60) 

This is our main analytic result. It shows that, regardless of rotation or p, the entropy 
of the neutral black brane can be recovered from the strongly coupled brane-antibrane 
field theory. It is also worth noting that this result cannot be obtained just by dimen- 
sional analysis. 

Alternatively, we can think of this result as constituting a highly nontrivial rela- 
tionship between the thermodynamical properties of near-extremal black branes and 
those of neutral black branes. 

Before moving on to develop some more physics of the brane-antibrane model, we 
may ask what our conclusion about the mass renormalization may do to other parts of 
our brane-antibrane-gas model. We know that 

-1 (g^) 



^ " dmsG ■ 

Therefore if, as we have shown, the entropies match but the masses are not equal, then 
the temperature is also affected in the same proportion as the mass 

(9-p) 

TpT = Th2^(^ . (62) 

(Note that this does not spoil any of our previous assumptions.) 

One might wonder how the unusual thermal properties of neutral black holes and 
black branes, such as negative heat capacity, could possibly be understood from an- 
alyzing this Dp-Dp field theory model, which is based on ordinary super- Yang-Mills 
systems in various dimensions. The explanation of this point in the context of the 
Dp- Dp model is simple and was given for the conformal cases in T. Upon inspection, 
one finds that there is a correlation between the energy in the open string gas and the 
contribution to the energy coming from the tension of the branes. Actually they are 
proportional to each other. This would give rise to the following interpretation for the 
negative specific heat. Namely, that since the energy in the gases is proportional to 
the tension energy, if we add more energy to the open string gases then we must create 
more Dp-Dp pairs to maintain the proportionality. This requirement to make more 
Dp-Dp pairs makes the system cool down, meaning that the more massive the black 
branes the colder they get. As we see, the moral of the story is that 'normal' SYM field 
theories on the worldvolume of the branes and antibranes behave 'abnormally', because 
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the total number of degrees of freedom controlled by N is not a constant; rather, it is 
given thermodynamically by the entropy maximization scheme. 

So it is natural and important to ask whether the same kind of correlations also 
hold in our particular system of interest, i.e. DO-DO, and in an even broader sense, for 
a general Dp-Dp system. 

It is easy to show the linear proportionality (|48jl between the energy in the gas and 
the energy contribution from the brane tension. Using data from the optimization of 
N 

iV~^, (63) 

Tp 

where mpT = 2NTp + mgas , gives 

"T-gas ~ NTp . (64) 

It is intriguing that exactly the same behaviour was observed in the case of non- 
dilatonic branes (D3, M2, M5). However, the way the result of 1_ was obtained, 
comparing energy density and pressure in field theory and supergravity, does not apply 
to our p = case of particular interest as there is no pressure on a + 1 dimensional 
world volume. 

It is important to note that in our entire analysis, the black objects under study are 
neutral - nonperturbatively nonextremal. It is quite pleasing that the brane-antibrane 
model precisely explains the entropy of the neutral black branes and black holes, up 
to renormalizations of the mass and angular momenta that we computed. 



3.4 Horizon size 

It is interesting to ask whether the transverse fluctuations of the Dp-Dp system in the 
microscopic picture can reproduce the size of the horizon of the corresponding black 
brane geometry. We now do a scaling analysis, not keeping precise numerical factors. 

Before addressing the general-p cases, it is instructive to review the simplest case 
p = 3. For near-extremal D3-branes at strong coupling, power counting in 
and conformal symmetry tell us that^ (X^)rms ~ NT"^. Using the model of to 
find the optimal value of gives ~ ttt-ft/'^s 

~ TTigas / T3 . Here we have used the 
information from the brane-antibrane model both to set the optimal value of A^ and 
to learn that there is roughly the same amount of energy density in the branes and 
in the open-string gas on those branes. Next, we use the assumption m-sc ~ ""t-fTj 
and the supergravity relationship between the Hawking temperature and the energy, 
'^SG ~ Tq/G. We also recall the assumption that the brane gas temperature is equal to 
the antibrane gas temperature, and both are equal to the Hawking temperature of the 
black brane. Collecting these facts, and using the relationship Gt| ~ 1, then gives a 
relationship between the supergravity horizon radius and variables in the field theory. 
The last assumption to be used is the supergravity relationship between the Hawking 
temperature and horizon radius Th ~ l/?'o- So {X'^)ms ~ ~ '''o- Therefore, 

the r.m.s. extent of the position fields of the field theory corresponds to the horizon 
radius. In other words, we started from a near-extremal gauge theory and ended up 
with the size of the horizon of a neutral black 3-brane. 

"^The r.m.s. expectation value is normalized with a factor of 
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Now we turn away from the conformal case. In later sections, we will be particularly 
interested in the p = case, so let us look at it here in some detail. 

In the mean field approximation, it has been shown JHj that the extent of the 
ground state of the DO-brane theory at low temperatures is controlled by t'Hooft cou- 
pling, not the temperature. In supergravity, on the other hand, the rough estimate 
of the size is temperature dependent. Therefore, the above logic that we used for the 
D3-brane case will not work for the DO-brane case. In |15j . a remedy for this confusing 
situation was proposed. Namely, that fast-varying degrees of freedom are physically 
inaccessible to a local supergravity observer. The supergravity probes just cannot re- 
solve high-frequency fluctuations of order of t'Hooft energies, that are occurring at 
a "microscopic" level in quantum mechanics. In other words, a supergravity probe 
cannot resolve distances sharper than /probe ~ /3. As a consequence, the supergrav- 
ity probes simply miss all the quantum dynamics of the DO-brane dynamics involving 
energies much higher than the Hawking temperature. Therefore, in making a micro- 
scopic model to reproduce supergravity, caution must be used regarding which degrees 
of freedom have to be included in the picture. 

Practically, the suggestion of ^Sj amounts to imposing a temperature-dependent 
cutoff on the spectral density of the transverse fluctuations, i.e. propagators, that 
leads to a temperature dependent {X'^)rms which qualitatively shows the same behavior 
as Tq does on supergravity side. 

So in what follows we will assume that, if one tries hard from the gauge theory 
point of view, the extension of the wavefunction of the near-extremal geometry should 
become consistent with supergravity expectations. We use this intuition motivated by 
IMSY in what follows. 

We now apply this thinking that came from studying DO-brane physics to the 
general-p case. We put together a number of ingredients that we have discussed in this 
section. We begin with 

mpT ~ ■m-branes + "^gas , (65) 

and note from the model of U we have 

'^gas ~ ^branes • 

(66) 

(Of course, the same energy density is in the brane gas and antibrane gas, so for the 
purposes of scaling we do not need to compute each contribution separately.) 

From e.g. [H], we have for the spatial extent of a near- extremal Dp or Dp geometry, 

Uo ~ (5^MA^)'/^'"^^r'/^'"^^ , (67) 

where Uq = r^/i'^. (Thinking of Uq as distance and T as energy, we see that this is 
exactly the same energy/distance relation found by ^1] for a D = 10 supergraviton 
probe of a Dp-brane geometry in the decoupling limit. The coincidence is not too 
surprising.) We now add in dynamical information from the brane-antibrane model, 

N ~ mpT/Tp , (68) 

to give 

r^-P ~ GmFTr2/(5-p) . (69) 
Lastly, we bring in the (fleld theory) relationship between temperature and mass 

T ~ (GmFT)"^/^^"^^ • (70) 
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Finally, we find 



(71) 



Since mpT ~ "msG , this is the same radius as that of the horizon in the neutral geometry 
of interest. Therefore, the brane-antibrane model does successfully provide a consistent 
picture of the spatial extent of the horizon. 

We now turn to checking the consistency of the supergravity approximation itself. 

3.5 Validity of supergravity 

The black Dp-branes on which we concentrate here are neutral. In order that they can 
be thought of as bona fide supergravity entities, we require at a minimum that the 
curvature at the horizon radius should be small in string units, to keep a' corrections 
small. In scaling, 

(Gm)i/(^-P) » £s . (72) 
Using the equilibrium value of N, this becomes 

(5s'^fA^Tp)i/(^-P)»4, (73) 

which in turn leads to 

<7,Af > 1 . (74) 
Of course, we also require that string loop corrections be under control 

9s <^ I, (75) 

for the D = 10 neutral black brane spacetime. 

Now, let us recall one important piece of physics from pP, the study of the conformal 
case. Even though the neutral geometry of interest is not studied in the decoupling 
limit, the brane and antibrane systems in the microscopic model are actually taken to 
be decoupled - as facts from AdS/CFT are used to describe the strongly coupled D3- 
brane theory. Therefore, at this point, it is appropriate to check whether the conditions 
for validity of the D = 10 near-extremal Dp-brane geometry are satisfied here also for 
the non-conformal C3fSGS, clS these geometries dictate for us the behaviour of the strongly 
coupled field theories on the Dp and Dp. 

Let us therefore study the IMSY conditions carefully, to gain understanding of the 
brane-antibrane side of the picture. As described in P] , there is a region in which type 
II D = 10 supergravity is valid 

1 « 5e'ff(f^) « , (76) 

where the effective coupling at energy scale U is 

gUu) ^ g'yuNUP-^ (77) 

Also, the horizon radius for the near-extremal geometry is related to the energy above 
extremality rrigas and the temperature T by 

Uo ~ (5YM"^gas)(^ ~ ig^y^N)(^T(^ . (78) 
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Now, our big neutral black brane will have a low Hawking temperature. This 
means that we are not in danger of violating the a' (left-hand) end of the bound ()76() . 
However, precisely because we are operating at such low temperatures, we may be 
concerned about violating the strong-coupling (right-hand) end of the bound (|76() . For 
example, for p = 0, we may be concerned about having to lift up to D = 11. 

It is a satisfying fact that our microscopic brane-antibrane model remains consistent 
in the Z) = 10 picture. The essential reason for this is that, in the brane-antibrane 
model, N is not an independent variable. Rather, N is thermodynamically determined. 
Using the fact that rrigas ~ NTp and formulae H77I78|) . and requiring that both ends of 
the IMSY bound (f7n|l are respected gives two conditions, which reduce to 

gsN->l, gs<.l- (79) 

As we saw before, open strings are strongly coupled, but closed strings are weakly 
coupled. 

Therefore, we do not need to be concerned about departing from the regime of 
validity of D = 10 supergravity, for our systems of strongly coupled branes and an- 
tibranes. 

In particular, for the p = case, we are always in the Z? = 10 supergravity regime. 
This means, in particular, that we never get to the D = 11 supergravity regime. 
Therefore, it is not apparent whether there is any simple relationship between these 
microscopic models for D = 10 neutral black branes and the BFSS matrix theory 
models, whose microphysical description is in terms of DO-brane degrees of freedom 
representing D = 11 M theory. 



4 Numerically investigating DO— DO physics 

For the remainder of this paper we will concentrate on the p = case, i.e. the brane- 
antibrane model of the = 10 Schwarzschild black hole. The reason is that the 
physics is a quantum mechanics, lending itself to the possibility of actually computing 
the behaviour of the strongly coupled field theory. 

4.1 Scalings, and strategy 

Of course, the idea of doing direct numerical simulations in the strongly coupled QM is 
not new. Kabat et al. [2] have used a method called Variational Perturbation Theory 
(VPT) to check the IMSY gravity/gauge duality conjectures [5] for the decoupling 
limit of Dp-brane systems with sixteen supercharges. In particular, the entropy of 

near extremal DO-branes was computed using VPT, and was found to match the 
with supergravity result to the level of approximation used [Hj. The agreement is 
impressive. In particular, the highly nontrivial temperature dependence of the free 
energy was obtained (approximately): f5F oc T^'*^! Indeed, this could be considered 
as the first (approximate) nonperturbative check of any gauge theory/gravity duality. 
The causal structure of spacetime from the point of view of the gauge theory also has 
been studied in the context of mean field Gaussian approximations . 

The next step for our program to understand = 10 Schwarzschild black holes 
would be to justify the crude brane-antibrane picture drawn in the last few sections, by 
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computing the microscopic entropy using the effective quantum description of a system 
of DO-DOs. 

Now, in the DO-DO system we cannot take a clean decouphng hmit, if we expect 
to keep the open string tachyon and massless string modes but not massive string 
modes. Therefore, our further progress in developing the brane-antibrane model is to 
be thought of as an approximate description, where the massive string modes are not 
fully decoupled. Of course, the same story was true of the Danielsson et al. work pp. 

Clearly, there are two different energy scales set by two different dimensionful cou- 
pling constants in our problem: a' and gi^yi = gs{27r)~'^ij^N . In gauge theory, every- 
thing is governed by ^y^, while a' controls the mass of the tachyon. For the physics 
of DO-branes alone, dimensional analysis and power counting in N tells us that a 
dimensionless quantity like PF can be written as 

Here, with the tachyon, we have two different scales: the t'Hooft coupling and a'. 
Thus, the free energy can be written as 

PF^o-m « ^'^^J^J^^'^'^aUNf') = N'g{^-^^^^, {gsNf/') . (81) 

Therefore, the free energy written in dimensionless units will be a function of QgN. 
Of course, it will also depend on the dimensionless inverse temperature measured in 't 
Hooft units, (3. 

In the limit where massive open string excitations can be neglected, and the string 
coupling is small, the system is relatively simple. We have two copies of the DO-brane 
theory - one for DO-branes and one DO-branes - plus a complex tachyon T and a 
massless Majorana fermion ^ coming from D-D open strings. 

The strategy is to compute the free energy of this system as a function of the 
tachyon classical background field, plus other quantities at finite inverse temperature 
p. In the field theory we compute at strong couplings using VPT, and in the closed 
string description the D = 10 supergravity approximation is valid. Then the following 
pieces of information could be read off immediately from the free energy 

Tachyon static mass The thermodynamically favourable value for the tachyon ex- 
pectation value can be computed by minimizing the effective action, i.e. the free 
energy of the system as a function of the tachyon background. 

Sign of the Tachyon dynamical mass Another important piece of information en- 
coded in the free energy is the effective mass of the tachyon which is given by 

where (T) is the tachyon expectation value. This effective mass includes the 
contributions from infinite numbers of loop diagrams via the Schwinger-Dyson 
equation. 

Magnitude of the Tachyon Dynamical Mass This quantity is a measure of the 
smallness of the tachyon fluctuations. A large dynamical mass results in small 
contributions to the entropy from the tachyon. 
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Phase Portrait of the theory The final goal is to calculate the phase portrait of 
the system, as a function of inverse temperature [5 and QsN . A sign change in 
m|, from negative to positive would be of great interest, since it would signal the 
tachyon stabilization phenomenon. This stabilization would give a justification 
for why a gas of DO and DO branes does not simply annihilate all the way down 
to a bunch of closed strings. 

One note of caution. The Hagedorn phenomenon in string theory might impact us 
here in one place: at temperatures of order of the string scale, massive open string ex- 
citations become important, but we are not incorporating massive excitations into our 
dynamics. In order to have a self-consistent description of the phenomena, therefore, 
we will look for a possible finite-temperature tachyon stabilization at a temperature 
well below the Hagedorn temperature. 

We now turn to describing the technology which we will use to perform the numer- 
ical simulation of the DO-DO system. 

4.2 Variational perturbation theory (VPT) 

The basic idea of VPT is simple. Using VPT to simulate strong-coupling DO-brane 
physics was first introduced in [21 El , and we review the salient points here for refer- 
nence. Suppose that the theory of interest has action 5, and the aim is to approximate 
its free energy at strong couplings. The idea is to use a free theory, with action Sq, with 
arbitrary tunable parameters. The next step is to find values for these parameters in 
such a way the free theory is a best fit for the full interacting theory. 
For any arbitrary Sq the following identity holds 

/?F = /?Fo - {e-^s~s,) _ 1^^^^ ^ (83) 

where refers to the free theory and C stands for connected contributions. Expanding 
the identity, 

(3F = (3Fo + {S- So)o -^iiS- Sof)o,c + ■■■ ■ (84) 

It is important to note that this expansion is not a perturbative expansion in the 
couplings of the original interacting theory of interest. 

Now, if terms to all orders in the above expansion are kept, then there is of course 
no dependence on the parameters of the trial free theory, as it is an identity true for 
any Sq. Taking a practical approach and terminating the series at any finite order, 
however, the series depends on the variational parameters. 

The next step is to fix the variational parameters. Minimizing the free energy with 
respect to those parameters (which is equivalent to requiring the trial free action to 
satisfy the Schwinger-Dyson equation) leads to a set of algebraic coupled equations 
(in general infinite in number), called "gap equations". Solution of the gap equations 
yields the variational parameters, which are then substituted back to obtain the free 
energy. 

This VPT method has in fact been checked explicitly for quantum mechanical sys- 
tems where calculations in the full interacting theory can actually be done exactly, and 
the above expansion captures the strong-coupling behaviour accurately; in particular, 
convergence is very fast. Our system of interest here is of course not solvable, so we 
need to use the approximate expansion procedure outlined above. 
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Practically speaking, of course, the infinite set of equations cannot be solved. There- 
fore it is necessary to cut off that infinite set of equations, and find the common roots 
of the finite coupled set of algebraic equations. Solving this system of equations is, at 
any rate, computationally very expensive. 

The last step is then to substitute the solution of the gap equations back into the 
free energy, to get the sought-after dependence on the parameters of interest. For us, 
these parameters will be the dimensionless inverse temperature /3, and the open-string 
coupling gsN. 

As outlined in [210], the VPT method is not straightforward when dealing with 
supersymmetry and gauge theory. Regarding gauge theory. Proposing a Gaussian 
theory for a non-dynamical gauge field in 0-1-1 dimensions is delicate, as one cannot 
simply gauge it away at finite temperature. Gauging away the gauge field leaves an 
observable (a Wilson loop), made from the zero mode of the gauge field encircling 
the Euclidean time direction, as a remnant. The Gaussian theory to be used is 
the one-plaquette model studied by Gross and Witten ^7j. Taylor-Slavnov identities, 
which are consequences of gauge symmetry and relate various correlation functions, 
get invalidated by the VPT expansion. Regarding supersymmetry. It turns out to 
be impossible in general to come up with a Gaussian theory with general variational 
parameters respecting supersymmetry, without inclusion of a trial action for auxiliary 
fields. A way to get around the problem is to work with the off-shell superspace 
formulation. 

On the other hand, the VPT method has some special positive features as well. 
One is that VPT automatically respects 't Hooft counting. Another is that VPT 
automatically cures infrared problems arising from having an infinite moduli space of 
vacua where DO-branes are far apart. The contribution coming from the zero mode 
sector goes like 0{N) and is therefore distinguishable from 0{N^) contributions. 

4.3 Motivating the action 

According to JSl (see also jl9j for a similar results), the action for a system of D9-D9 
with the U{1) gauge fields living on its world- volume is given by the following^ 



where we have used 'bars' to denote quantities in the DO sector and 'no bars' to denote 
quantities in the DO sector. Of course, T is a complex bi-fundamental gauge field so 
that Dfj^T = d^T — iTA^ + iA^T. Dimensional reduction of this action should give 
the correct action for the dynamics of lower dimensional Dp- Dp systems including our 
p = case. As usual, we substitute Ai = (27ra')~^Xj , Ai = (27ra')~^Xj. Performing 

^We are using the mostly plus convention for the metric signature. 



S 



2Tg / d^^x exp(-27ra'rtr) 1 + Svra'^ ln(2)L>^rtL>^'r + 




(85) 
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field redefinitions, a nonabelian generalization of this toy approximate action is 



S = — dx^Tv exp 



41n2 J 



^ {doT^d^T + idoT\A°T -TA^) 



+i(A°T"f - T^A^)doT - (A°rt - T^A^){AoT - TAq 
-(X^rt - T^T)(X^T - TXi)) + 2roIiVx^ 



+5U[^^^°] + 5W^,^Y (86) 

where we use the usual definition of the 't Hooft coupling, I/^ym ~ ^Tr^gi"^^^, and 
where S'qq and are the corresponding low-energy actions for DO and DO respec- 
tively. Obviously, T is in the N x N bifundamental representation. 

Now, in order to be able to compute the thermal partition function, we need to 
Euclideanize the Minkowskian action. First, we analytically continue the timelike 
component of both U{N) gauge fields as well as Minkowskian time direction while 
leaving the other fields untouched, 

xO = -ir, iSE = SM, iA% = A%, z = ^4^ • (87) 

Next, we Fourier expand T, A^, A^ , X* and X* in r with a periodicity /? given by 
the the inverse temperature. Since all these variables are bosonic, they have periodic 
boundary conditions, 



AoiT),Ao{T),X\T),X\T),T{T)] = ^ E {^oi, Aoi, Xi Xl T,} 



Writing the action in terms of Fourier modes is straightforward. We take the above 
expansions and perform the Euclidean time integral, using the following variational 
correlators in Wick contractions of Feynman diagrams contributing to the free energy 
(in our case, the amplitude {She — Soe)o,) 

{^eAB^LcD)o = ^al^"'''^e+mSADSBC a,6 = 3..9, 
(^OOAB^OOCd)o = Po^AD^BC , 

{'^eAB'^mCD)o = iJ^lm^AD^BC ■ (89) 

and similarly for the barred variables. The need for separation of transverse scalar 
directions comes from the limitations of the original superspace formulation [Hj of the 
DO-brane super symmetric matrix quantum mechanics. Namely, for each copy of the 
QM, there are 2 scalars X* in the gauge multiplet and 7 other (j)"" scalars coming 
from scalar multiplet. In other words, the original 50(9) i?-symmetry is broken to 
50(2) X Ga {(t)" are a 7 of G2). 

The numerical simulation of this action, which is an approximation to the system of 
DO-branes, anti-DO-branes and the tachyon sector where both copies of the DO-brane 
theory are supersymmetric is computationally very expensive, partly because of the 
need to find the solution to hundreds of nonlinear equations. In addition, we have to 
cover mapping out of the phase portrait as a function of (3 and QsN . We therefore 
report on some initial results in a more stripped-down version of our model. 
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4.4 The toy model 

In this section, we look at a "toy" model which consists of large- bosonic matrix 
quantum mechanics (representing the DO-brane and anti-DO-brane theories) plus our 
bi-fundamental charged tachyon. 

This toy model is obviously a good approximation to the original problem in the 
high-temperature limit where the Euclidean time circle is so small that nonzero thermal 
KK modes are heavy. In this limit, there is no antiperiodic fermion left after reduction 
to -|- dimensions. Of course, we are interested in low temperatures where the 
supergravity description is valid and quantum mechanics is strongly coupled. 

There are reasons to believe why this toy model could demonstrate some generic 
behaviors of interest to us, which will continue to be true even in the case where we have 
two copies of supersymmetric matrix quantum mechanics. It is well known |^ that the 
extent of the ground state of a system of DO-branes, defined in terms of fluctuations of 
transverse scalars, is controlled by the 't Hooft coupling. It is intuitively plausible to 
expect that, at strong coupling, the tachyon could acquire a large induced mass through 
its coupling to these large transverse scalar fluctuations. As has been observed in |2j, 
both bosonic matrix quantum mechanics and its supersymmetric version have large 
transverse fluctuations. So if one takes this intuition seriously, as far as the dynamics 
of the tachyon is concerned in the low energy approximation, even the bosonic version 
could be a good approximation to the right physics^ 

In the model of Danielsson et al (1^ which motivated our work here, the dynamical 
assumption is that the tachyon field is stabilized by finite-temperature strong-coupling 
effects about T = 0. We therefore expand T = + t, where t is the fluctuation part, 
and plan to show the consistency of this assumption a posteriori by showing that the 
mass is indeed positive - and large^ - at strong coupling. Therefore, the terms relevant 
to our discussion which are quadratic in the tachyon fluctuations t can be written 



Se = I dT{\TTDrX'D^X' -\t,[X\X^][X\X^]) 

aW Jo 2 4 

+ -r- I dTi^TlDrltDrT -\TT\X\X'^]\jt,X^]) 

9ym Jo 4 

+ [ dTTr[^— (drt^drt + A'^thA'^ + tt:4°:4°t + xhUx' + t^TTt 

Jo ^^YM V 



(90) 



Morally, we are going to treat the tachyon as background. By this we simply mean 
that we are interested in finding the free energy as a function [al] of the tachyon, after 
integrating out the other fields in the problem - nonperturbatively using VPT. For this 
reason, no trial action for the tachyon is introduced; the VPT techniques are used for 
the fields other than the tachyon. 



^Although irrelevant to the above discussion, the difference is in the actual value of the ground state 
energy, which is - for each copy - of course zero in the supersymmetric case but large and positive in the 
bosonic case. 

^We have ignored the higher powers in the expansion of the exponential prefactor in (|86(l sitting in front 
of the Lagrangian; this will be justified a posteriori by the large-mass finding. 
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Using VPT, we have 

/3F = (3Fo + PFo + (5do + Sj^ + St- So)o - ^((5m)')o,c (+•..), (91) 

where St is the tachyon contribution to the action, and the subscript zero refers to the 
contribution coming from the trial free action for the other fields. 

As a first step, we can neglect the terms in the action that are cubic in the tachyon, 
as far as the VPT expansion is concerned. There is a good reason why this is justified. 
Namely, because the cubic contribution to the free energy enters at the order of 0{t'^), 
which is small at small-t. 

Following the story for a single clump of DO-branes pj| , we put in the ghost field from 
gauge-fixing, whose Fourier modes are denoted s^^o (obviously, for physical reasons 
there is no zero mode) . We also use equations from pages 25-27 of that work; to save 
space here we only write here what changes when we introduce our tachyon field. 

We also introduce convenient dimensionless units. These are defined such that in 
term of them one gets an overall factor of A^^ in the free energy, and the rest is just a 
function of dimensionless quantities and gsN. We define newly dimensionless variables 
with tildes; we have 

and similarly for the barred variables. Here / = (gyj^N)^/'^ . Note that in this bosonic 
model we do not have to treat the X*, X* fields for i = 1, 2 differently than for i = 3 . . .9 
because there is no need to worry about maintaining explicit super symmetry; therefore 
we do not have and we just use to represent correlators of all nine position fields. 
Note that we are measuring the tachyon expectation value in string units. 

Henceforth, we drop the tildes for notational simplicity. Then we obtain the fol- 
lowing expression for the free energy as a function of the tachyon 

PF = /3F(A)n + /3F(A)n + ^(Tr(C/ + C/t))n + ^(Tr(I7 + I7^))n 
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^ log logaj + N^}_^ log se + N^}_^ log 

9iV2^ 27r/2 2 ..,9A^2 2^, 

E((^)'^^ + 1) + E((^)'^^ + 1) + E 
+^(E^^ + E'"^)ETr 



2/3 ^Z^"'^ ' A^-^'A^-N^ 



81n2(5,iV)2/3 Ar2 {g.Nfl^ 



22 



In what follows, we use a natural ansatz for the propagators 



-l=[{^-f? + m') , (94) 

and similarly for the barred variables. 

The gap equations are easily obtained. We have (and similarly for barred variables) 

where po can be explicitly obtained in terms of A, the coupling in the one-plaquette 
action, because the one-plaquette theory is soluble in the large- A/" limit. 
For A < 2, 




(l-^)log(l-^) 



0, (96) 



while for A > 2, 



We will solve this system of highly nonlinear, coupled equations. The strategy will 
be to treat the tachyon as a background. Solving the gap equations, and substituting 
back the solution as a function of the tachyon, then gives the free energy as a function 
of the (background) tachyon. Given this effective action for the tachyon, it is then 
straightforward to read off the tachyon effective mass. The plan will be to work out 
the two-dimensional "phase portrait" of the system (with coordinates ggN and (3), and 
find out whether any tachyon stabilization is possible. By this we simply mean that 
there are points in this two-dimensional portrait where the effective tachyon mass is 
positive. 

Before we move to displaying our numerical results, a number of remarks are in 
order. 

• As we saw earlier, in the supcrgravity limit in which wc arc interested, QsN S> 1. 
By inspection of the equations it might then seem that, in the units we use, 
the wrong-sign mass term for the tachyon is suppressed and we are done! This 
quick-and-dirty reasoning is, however, too naive. In fact, in order not to excite 
any massive open string state, the temperature must be small in string units. In 
turn, this implies that 

/?»(<7sA^) (98) 

where /? is the dimensionless inverse temperature. With such low dimensionless 
temperatures, it becomes more and more difficult to overcome the wrong-sign 
mass term through the tachyon couplings to the transverse coordinates 
and the world- volume gauge fields ^o- 

• It is crucial to remember that all variational parameters are implicit functions of 
the tachyon t. So it is certainly not true that we can simply read off the effective 
tachyon mass from the the coefficients of the quadratic tachyon fluctuation in the 
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expression for the free energy. For example, terms like {N'^ /P)Po'^i(^i, which 
do not have any explicit background tachyon dependence, are becoming tachyon- 
dependent through the gap equations. All such terms make contributions to 
the effective tachyon mass, as one expands the variational parameters around 
^£ Trtjt^ = 0. (Note that these expansions exist, since the solution to the gap 
equations is differentiable with respect to t in the vicinity of ^^Trt^t^ = 0.) 

• Obviously, it is always possible to stabilize the tachyon at weak couplings, but 
the temperatures needed for this are above the Hagedorn temperature. It is easy 
to see this in perturbation theory, partly because at such high temperatures the 
DO QM is perturbative. Perturbatively, to stabilize the tachyon, it is sufficient to 
have 

for small gsN. 

In our investigations, we use the Metropolis algorithm. This is a fast method gen- 
erally; it really comes into its own for numerical investigations of the supersymmetric 
model, the numerical results from which we hope to report on in the near future. We 
summarize some important features of the Metropolis algorithm in the Appendix. 



4.5 Numerical results 

We are interested in the sign of the mass of the tachyon zero mode, for each point 
in the two-dimensional parameter space, coordinatized by QgN and /? which is the 
dimensionless inverse temperature in 't Hooft units. 

We now plot our results for the phase portrait of the toy model theory. Inverse 
temperature, in dimensionless 't Hooft units, is on the y-axis, while the open-string 
coupling gsN on the x-axis. (As a reminder, the parameter QgN tells us the bare 
negative mass-squared of the tachyon.) Points at which the tachyon mass squared is 
positive are indicated by a cross, and those at which it is negative are indicated by 
a circle. The green curve delineates where the Hagedorn phenomenon is important; 
points above the curve are below the Hagedorn temperature. The red straight line 
delineates the region where the gauged QM theory has a strong effective coupling 
(set by the inverse temperature and 't Hooft coupling); the temperatures where the 
quantum mechanical theory is strongly coupled lie above the red line. 

For clarity, we separate the regimes of smaller ggN from the regime of larger ^^A^. 

As we expected for tachyon stabilization, at small QgN, one has to go to very high 
temperatures; this is clear from Figure ^ As we see for large QsN, on the other hand, 
even for low temperatures compared to the Hagedorn temperature, there is a possibility 
of tachyon stabilization. 



5 Discussion 

In this paper, following the lines of Jl, we have examined the validity of the brane- 
antibrane model for neutral black branes of various dimensionality. What we have 
found is that this Dp-Dp picture works for arbitrary p. Interestingly, this is true 
even for the "dilatonic" branes, where the supergravity entropy cannot be written in 
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Figure 1: Plot of results for the ggN < 2 section of the phase portrait. 




Figure 2: Plot of results for the ggN > 2 section of the phase portrait. 
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terms of the entropy of a free world- volume gauge theory. In the non-conformal cases, 
using IMSY duality proves to be working, relating strongly coupled gauge theories to 
the corresponding near-extremal backgrounds. We were able to show that, even after 
adding general angular momenta quantum numbers, the brane-antibrane model is able 
to reproduce the supergravity entropy. The exact agreement in the entropy between 
the microscopic brane-antibrane side and the supergravity side comes at the expense 
of a renormalization of the mass by a p-dependent power of two, which we interpret as 
a binding energy, and a renormalization of the angular momenta by a simple factor of 
two. We also find that the extent of the wavefunction of the Dp — Dp system exhibits 
the same scaling with the mass as does the radius of the horizon of the corresponding 
black brane. 

In the original proposal of the authors argued that the tachyon stabilization 
at finite temperature and strong coupling could lead to a reappearance of open string 
degrees of freedom. Even at the field theory level, questions regarding strong coupling 
dynamics are often hard to answer, but in the special case of DO-DO there is a possibility 
to address some of the aspects in the context of a simple toy model following the 
lines of [210] • Using techniques developed there, we found the effective action for the 
tachyon numerically, in the context of the toy model. We have seen signals of tachyon 
stabilization in this toy model, at strong open string couplings and low temperatures 
compared to the string scale. This might open a window to see clearly why it is possible 
to have a long-lived state of DO-branes and DO-branes, without decaying into closed 
strings at couplings of order one and low temperatures. 

Our investigation here was for the case of neutral black branes, where the numbers 
of microscopic constituent branes and antibranes A'^ are equal. A natural extension 
of our work here will be to consider the case where the D = 10 black branes carry 
p-brane charge Qp = N — N. We hope to report on this story in the future®. 

Inspired by numerical investigations, we might speculate about the possible shape of 
the potential^ in this model, say for the simplest case of a diagonal tachyon condensate 
at finite temperatures 

PF oc e-^l'^l'(l + X{N,gYM,f3)\T\^ + . . .) (100) 

where A > 7. This potential has one minimum and two symmetrically located maxima. 
The minimum would correspond to the tachyon stabilization, i.e. the tachyon particle 
is classically stable^", even though this vacuum would be unstable due to quantum 
and thermal fluctuations and the finiteness of the height of the barrier. The tachyon 
can be expected to decay eventually into closed string radiation. In other words, 
this temporary and approximate decoupling between the open string modes and the 
closed strings degrees of freedom - which has apparently manifested itself in our ability 
to compute the entropy of the closed-string background out of Yang-Mills degrees of 
freedom in the context of brane-antibrane model - would not last forever, and the 
tachyon particle ground state would have a finite lifetime. If this picture corresponds 
well to the true stringy physics, then it would be plausible to think of the unstable 

^The p = 3 case with charge was of course previously considered in [n i21j . It will be interesting to 
whether this agreement holds up also in the noM-conformal cases. We thank Alberto Guijosa for a disc 
in this regard. 

^at least for |r| in the vicinity on T = 

"'^'^The kinetic term probably is not canonically normalized in the original field. 
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maxima discussed above as being the unstable state of a DO-DO pair created from the 
energy in the gas Uving on the world volume of the clump of DOs and DOs, due to the 
quantum and thermal fluctuations. This interpretation also leads to an estimate for 
the lifetime. Climbing up the hill to create such a pair requires energies at least of the 
order of 1/ gs, so the whole tunnelling process is suppressed by e"^/^". 
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Appendix 

The problem of finding the common solution to a set of nonlinear coupled algebraic 
equation can be though of a minimization process. To see this, consider the collection 
of the equations to be solved, Q = {/i = 0, ■■.,fn = 0}. Clearly, a quadratic form 7i 
can be made out of this, viz. Ti = Ylifi- there is a solution to G, this solution 
would be the global minimum of Ti. i.e.,/i = 0, /„ = 0. 

To solve this nonlinear system of equations, we use the Monte Carlo local update 
method ( "Metropolis" algorithm) . This method is one of the stochastic search methods 
widely used for minimization of multivariable functions. 

This method does not get trapped in local minima, because of a smart choice of 
control variables. The way this method accomplishes the goal is that, by introducing an 
unphysical temperature, thermal fluctuations of the configuration vector are allowed 
to happen in all possible directions in the configuration space - even in the "wrong 
directions" (which increase the Hamiltonian) . By decreasing the temperature very 
slowly, one lets the vector find the global minimum (or minima). For a vector trapped 
in a local minimum, there would be a certain probability to climb up the barrier 
because of the thermal fluctuations, and so the vector could find its way all the way 
down to the global minimum subject to a large number of iterations. This story is 
actually very similar to the physics of the annealing process, and it is sometimes called 
the "stimulated annealing" method. With this method, a vector stores the starting 
point configuration. This configuration is changed randomly to generate new child 
configurations. A new configuration gets accepted with the probability 1 if it has 
a smaller Hamiltonian than its parent, but if it has not then it would be accepted 
with the following Boltzmann-like probability -Pacceptance ~ Q-l^fakc^T-i^ where /3fake is 
a unphysical parameter acting like an inverse temperature, and ATC is the energy 
difference between the child and parent configurations. After many iterations, and at 
the same time lowering the temperature slowly, the configuration vector stabilizes on 
the one(s) which minimize(s) the Ti.. 

The Metropolis algorithm has strong advantages over the Newton-Raphson (NR) 
method, especially when the dimensionality of the configuration space is big. Under 
such circumstances, the NR method becomes very inefficient, and its starting point 
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dependence can be prohibitively large. Practically, this makes it very hard to reach the 
solution in reasonable computer time. The great advantage of the Metropolis method 
is the fact that it is almost independent of starting point. Another advantage of 
Metropolis is that it naturally avoids getting stuck in local minima, a problem to 
which NR often succumbs. 
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